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Abstract. We consider a general theory of curvatures of discrete surfaces equipped 
with edgewise parallel Gauss images, and where mean and Gaussian curvatures of 
faces are derived from the faces' areas and mixed areas. Remarkably these notions 
are capable of unifying notable previously defined classes of surfaces, such as discrete 
isothermic minimal surfaces and surfaces of constant mean curvature. We discuss 
various types of natural Gauss images, the existence of principal curvatures, constant 
curvature surfaces, Christoffel duality, Koenigs nets, contact element nets, s-isothermic 
nets, and interesting special cases such as discrete Delaunay surfaces derived from 
elliptic billiards. 



1. Introduction 

A new field of discrete differential geometry is presently emerging on the border 
between differential and discrete geometry; see, for instance, the recent books (2], [6]. 
Whereas classical differential geometry investigates smooth geometric shapes (such as 
surfaces), and discrete geometry studies geometric shapes with a finite number of el- 
ements (such as polyhedra), discrete differential geometry aims at the development of 
discrete equivalents of notions and methods of smooth surface theory. The latter appears 
as a limit of refinement of the discretization. Current progress in this field is to a large 
extent stimulated by its relevance for applications in computer graphics, visualization 
and architectural design. 

Curvature is a central notion of classical differential geometry, and various discrete 
analogues of curvatures of surfaces have been studied. A well known discrete analogue 
of the Gaussian curvature for general polyhedral surfaces is the angle defect at a vertex. 
One of the most natural discretizations of the mean curvature of simplicial surfaces (tri- 
angular meshes) introduced in |13] is based on a discretization of the Laplace-Beltrami 
operator (cotangent formula). 

Discrete surfaces with quadrilateral faces can be treated as discrete parametrized 
surfaces. There is a part of classical differential geometry dealing with parametrized 
surfaces, which goes back to Darboux, Bianchi, Eisenhart and others. Nowadays one 
associates this part of differential geometry with the theory of integrable systems; see 
[9], 117] . Recent progress in discrete differential geometry has led not only to the dis- 
cretization of a large body of classical results, but also, somewhat unexpectedly, to a 
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better understanding of some fundamental structures at the very basis of the classical 
differential geometry and of the theory of integrable systems; see [6]. 

This point of view allows one to introduce natural classes of surfaces with constant 
curvatures by discretizing some of their characteristic properties, closely related to their 
descriptions as integrable systems. In particular, the discrete surfaces with constant 
negative Gaussian curvature of |18] and |24] are discrete Chebyshev nets with planar 
vertex stars. The discrete minimal surfaces of [3^ are circular nets Christoffel dual to 
discrete isothermic nets in a two-sphere. The discrete constant mean curvature surfaces 
of |4j and [10 are isothermic circular nets with their Christoffel dual at constant distance. 
The discrete minimal surfaces of Koebe type in [TJ are Christoffel duals of their Gauss 
images which are Koebe polyhedra. Although the classical theory of the corresponding 
smooth surfaces is based on the notion of a curvature, its discrete counterpart was 
missing until recently. 

One can introduce curvatures of surfaces through the classical Steiner formula. Let us 
consider an infinitesimal neighborhood of a surface m with the Gauss map s (contained 
in the unit sphere S 2 ). For sufficiently small t the formula 

TTi = m + ts 

defines smooth surfaces parallel to m. The infinitesimal area of the parallel surface m t 
turns out to be a quadratic polynomial of t and is described by the Steiner formula 

(1) dA(m*) = (1 -2Ht + Kt 2 )dA(m), 

Here dA is the infinitesimal area of the corresponding surface and H and K are the 
mean and the Gaussian curvatures of the surface m, respectively. In the framework of 
relative differential geometry this definition was generalized to the case of the Gauss map 
s contained in a general convex surface. 

A discrete version of this construction is of central importance for this paper. It 
relies on an edgewise parallel pair m, s of polyhedral surfaces. It was first applied in 
|19|, 120] to introduce curvatures of circular surfaces with respect to arbitrary Gauss 
maps s £ S 2 . We view s as the Gauss image of m and do not require it to lie in S 2 , 
i.e., our generalization is in the spirit of relative differential geometry [22.1. Given such 
a pair, one has a one-parameter family m = m + ts of polyhedral surfaces with parallel 
edges, where linear combinations are understood vertex- wise. 

We have found an unexpected connection of the curvature theory to the theory of 
mixed volumes |21j . Curvatures of a pair (m,s) derived from the Steiner formula are 
given in terms of the areas A{m) and A(s) of the faces of m and s, and of their mixed 
area A(rn, s): 

A(m t ) = (l-2Ht + Kt 2 )A(m), H = K = 4r\- 

The mixed area can be treated as a scalar product in the space of polygons with parallel 
edges. The orthogonality condition with respect to this scalar product A(m, s) = 
naturally recovers the Christoffel dualities of [3J and [TJ, and discrete Koenigs nets (see 
[6]). It is remarkable that the aforementioned definitions of various classes of discrete 
surfaces with constant curvatures follow as special instances of a more general concept 
of the curvature discussed in this paper. 
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It is worth to mention that the curvature theory presented in this paper originated 
in the context of multilayer constructions in architecture |15] . 

2. Discrete surfaces and their Gauss images 

This section sets up the basic definitions and our notation. It is convenient to use 
notation which keeps the abstract combinatorics of discrete surfaces separate from the 
actual locations of vertices. We consider a 2-dimensional cell complex (V, E, F) which 
we refer to as mesh combinatorics. Any mapping m : i £ V \—> mi 6 M 3 of the vertices 
to Euclidean space is called a mesh. If all vertices belonging to a face are mapped to 
co-planar points, we would like to call the mesh a polyhedral surface. If / = . . . ,i n ) is 
a face with vertices i±, . . . , i n , we use the symbol m(f) to denote the n-gon , . . . , rrii n . 

Definition 1. Meshes m,m! having combinatorics (V,E,F) are parallel, if for each 
edge £ E, vectors nii — rrij and m\ — m!- are linearly dependent. 

Obviously for any given combinatorics there is a vector space (M 3 )^ of meshes, and 
for each mesh there is a vector space of meshes parallel to m. If no zero edges with 
rrii = rrij are present, parallelity is an equivalence relation. In case m is a polyhedral 
surface without zero edges and m! is parallel to m, then also m! is a polyhedral surface, 
such that corresponding faces of m and m! lie in parallel planes. 

A pair of parallel meshes m, m! where corresponding vertices nn, m\ do not coincide 
defines a system of lines Li = rrii V m^. By parallelity, lines associated with adjacent 
vertices are co- planar, so the lines Li constitute a line congruence [6]. It is easy to 
see that for simply connected combinatorics we can uniquely construct m! from this 
congruence and a single seed vertex m' io S Li , provided no faces degenerate and the 
lines Li intersect adjacent faces transversely. 

A special case of this construction is a parallel pair m, ml of polyhedral surfaces 
which are offsets at constant distance d of each other, in which case the lines Lj are 
considered as surface normals. The vectors 

Si = - mi) 

define the mesh s called the Gauss image of m. Following j!5l 116] . we list the three 
main definitions, or rather clarifications, of the otherwise rather vague notion of offset: 

* Vertex offsets: the parallel mesh pair m, m' is a vertex offset pair, if for each vertex 
i G V, || mj — ?7^|| = d. The Gauss image s is inscribed in the unit sphere S 2 . 

* Edge offsets: (m, m') is an edge offset pair, if corresponding edges mirnj and m\m^ 
are contained in parallel lines of distance d. The Gauss image s in midscribed to the 
unit sphere (i.e., edges of s are tangent to S 2 and s is a Koebe polyhedron, see [T]). 

* Face offsets: (m, m') is an face offset pair, if for each face f £ F, the n-gons m(f), 
m'(f) lie in parallel planes of distance d. The Gauss image s is circumscribed to S 2 . 

The polyhedral surfaces which possess face offsets are the conical meshes, where 
for each vertex the adjacent faces are tangent to a right circular cone. The polyhedral 
surfaces with quadrilateral faces which possess vertex offsets are the circular surfaces, 
i.e. their faces are inscribed in circles. 
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Remark 1. Meshes which possess face offsets or edge offsets can be seen as entities of 
Laguerre geometry |14] , while meshes with regular grid combinatorics which have vertex 
offsets or face offsets are entities of Lie sphere geometry [5j E] . 

3. Areas and mixed areas of polygons 

As a preparation for the investigation of curvatures we study the area of ra-gons in 
M? . We view the area as a quadratic form and consider the associated symmetric bilinear 
form. The latter is closely related to the well known mixed area of convex geometry. 

3.1. Mixed area of polygons. The oriented area of an n-gon P = (po, . . . ,p n _i) 
contained in a two-dimensional vector space U is given by Leibniz' sector formula: 

(2) A ( p ) = lj2 n< ^ det(pi,p i+1 ). 

Here and in the following indices in such sums are taken modulo n. The symbol det means 
a determinant form in U. Apparently A{P) is a quadratic form in the vector space U n , 
whose associated symmetric bilinear form is also denoted by the symbol A(P, Q): 

(3) A(XP + fiQ) = X 2 A{P) + 2\fiA(P, Q) + fi 2 A{Q) . 

Note that in Equation ^ the sum of polygons is defined vertex-wise, and that A(P, Q) 
does not, in general, equal the well known mixed area functional. For a special class of 
polygons important in this paper, however, we have that equality. 

Definition 2. We call two n-gons P,Q £E U n parallel if their corresponding edges are 
parallel. 

Lemma 3. If parallel n-gons P,Q represent the positively oriented boundary cycles of 
convex polygons K,L, then ([3| computes the mixed area of K,L. 

Proof. For A, \x > 0, the polygon XP + \iQ is the boundary of the domain XK + p,L, 
and so ([3]) immediately shows the identity of A(P, Q) with the mixed area of K, L. □ 

In view of Lemma [3j we use the name mixed area for the symbol U A(P, Q)" in 
case polygons P, Q are parallel. Next, we consider the concatenation of polygons P\,P2 
which share a common sequence of boundary edges with opposite orientations which 
cancel upon concatenation. Successive concatenation of polygons P%, . . . ,P& is denoted 
by Pi © ... © Pk- It is obvious that ^4(©j P-i) = ^A(Pi), but also the oriented mixed 
areas of concatenations have a nice additivity property: 

Lemma 4. Assume that Pi©- • -ffiPfc and P[(B - ■ '®Pi are two combinatorially equivalent 
concatenations of polygons, and that for i = 1, . . . , k, polygons Pi, P[ are parallel. Then 

(4) A{P X @ ■ ■ ■ © P k , P{ © • • • © P' k ) = A{P U P[) + ■■■ + A(P k , P' k ). 

PROOF. It is sufficient to consider the case k = 2. We compute A(P\ © P 2 , P{ © P2) = 
i\ t=0 A((P 1 ®P 2 ) + t(P[(BP^) = f t \ t=0 A((P 1 + tP[)®(P 2 + tPi)) = f t \ t=o A(P 1 +tP[) + 
I j t=0 ^(P 2 + tP> 2 ) = A(P U P{) + A{P> 2 , P' 2 ). □ 
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Figure 1: (a) Parallel quadrilaterals whose vertices lie on the boundary of their convex hull, 
(b) Parallel quadrilaterals whose vertices do not lie on the boundary of their convex hull. 



3.2. Signature of the area form. We still collect properties of the mixed area. 
This section is devoted to the zeros of the function A{xP -\- yQ), where P, Q are parallel 
n-gons in a 2-dimensional vector space U. 

Theorem 5. Consider a quadrilateral P which is nondegenerate, i.e., three consecutive 
vertices are never collinear. Then the area form in the space of quadrilaterals parallel 
to P is indefinite if and only if all vertices po, ■ ■ ■ ,P3 are extremal points of their convex 
hull. If P degenerates into a triangle, then the area form is semidefinite. 

Proof. We choose an affine coordinate system such that P has vertices (^), ([]), (J), 
(^) (cf. Figure flj). Translations have no influence on the area, so we restrict ourselves 
to computing the area of Q parallel to P with q x = Q, q 3 = (*). Then A(Q) = 
(s' t') ■ ^~^ s , • The determinant of the form's matrix equals (1 — s — t)/st, 

so the form is indefinite if and only if two or none of s,t, 1 — s — t are negative, i.e., all 
vertices lie on the boundary of the convex hull. In the degenerate case of three collinear 
vertices we compute areas of triangles all of which have the same orientation. □ 

Proposition 6. Assume that n-gons P, Q are parallel but not related by a similarity 
transform. Consider the quadratic polynomial <p(x,y) = A{xP + yQ). 

1. Suppose there is some combination P 1 = XP + uQ which is the vertex cycle of a 
strictly convex polygon K. Then <p factorizes and is not a square in M[x,y\. 

2. Assume that n = 4 and that some combination XP + uQ is nondegenerate. Then 
ip is no square inW[x, y]. It factorizes <^=^ the vertices of XP + uQ are extremal points 
of their convex hull. 

Proof. 1. Change (X,u) slightly to (A',//), such that |**'| / and Q' := X'P + u'Q 
still bounds a strictly convex polygon, denoted by L. Consider (p'{x, y) = A{xP' + 
yQ'). As (p and ip' are related by a linear substitution of parameters, it is sufficient to 
study the factors of ip': According to (|3]), the discriminant of ip' equals 4(A(P',Q f ) 2 — 
A(Q')A(P')) = 4(A(K,L) 2 - A(K)A(L)), which is positive by Minkowski's inequality 
|21j . The statement follows. 

In case 2 we observe that any element polygon parallel to P arises from some xP+yQ 
by a translation which does not change areas. It is therefore sufficient to consider the 
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areas of the special quads treated in the proof of Theorem |5] The matrix of the area 
form which occurs there is denoted by G. Obviously tp factorizes <^=^ det G < 
the area form is indefinite or rank deficient. We see that det G / 0, so rank deficiency 
does not occur (and consequently <p is no square). We use Theorem [5] to conclude that 
tp factorizes -4=>- the vertices of XP + [iQ lie on the boundary of their convex hull. □ 



4. Curvatures of a parallel mesh pair 

Our construction of curvatures for discrete surfaces is similar to the curvatures de- 
fined in relative differential geometry [22 , which are derived from a field of 'arbitrary' 
normal vectors. If the normal vectors employed are the usual Euclidean ones, then the 
curvatures, too, are the usual Euclidean curvatures. 

A definition of curvatures which is transferable from the smooth to the discrete 
setting is the one via the change in surface area when we traverse a 1-parameter family 
of offset surfaces. Below we first review the smooth case, and afterwards proceed to 
discrete surfaces. 



4.1. Review of relative curvatures for smooth surfaces. Consider a smooth 
2-dimensional surface M in R 3 which is equipped with a distinguished "unit" normal 
vector field n : M — > M 3 . It is required that for each tangent vector v 6 T p M, the vector 
dn p {v) is parallel to the tangent plane T p M, so we may define a Weingarten mapping 
<7p : T p M — * T p M by o~ p (v) = —dn p {v) (a unit normal vector field in Euclidean space 
M 3 fulfills this property). Then Gaussian curvature K and mean curvature H of the 
submanifold M with respect to the normal vector field n are defined as coefficients of 
cip's characteristic polynomial 

(5) Xa p (A,/x) := det(Aid + fia p ) = A 2 + 2XfiH(p) + fi 2 K{p). 

We consider an offset surface M s , which is the image of M under the offsetting map 
e 5 : p i— ► p + 8 ■ n(p). Clearly, tangent spaces in corresponding points of M and M s are 
parallel, and corresponding surface area elements are related by 

dA s 

(6) — = det(de^) = det(id + 6 ■ dn) = det(id - 5a p ) = 1 - 25H + 5 2 K, 
dA p 

provided this ratio is positive. This equation has a direct analogue in the discrete case, 
which allows us to define curvatures for discrete surfaces. 



4.2. Curvatures in the discrete category. Let m be a polyhedral surface with 
a parallel mesh s. We would like to think of s as the Gauss image of m, but so far s is 
arbitrary. The meshes m 5 are offsets of m at distance 5 (constructed w.r.t. to the Gauss 
image mesh s). For each face / € F, the n-gons m(f), s(f), and m s {f) lie in planes 
parallel to some two-dimensional subspace Uf. The area form in Uf and the derived 
mixed area are both denoted by the symbol A. We have the following property: 
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Theorem 7. If m, s is a parallel mesh pair, then the area A(m s (f)) of a face f of an 
offset m 5 = m + 5s obeys the law 

(7) A(m s (f)) = (1 - 25H f + 5 2 K f ) A{m{f)) , where 

( R ) H A(m(f),s(f)) = A(s(f)) 

{ ' f A(m(f)) ' / A{m{f)Y 

Proof. Equation ([7]) can be shown face- wise and is then a direct consequence of (|3]). 
As all determinant forms in a vector space are multiples of each other, neither Hf nor 
Kf depend on the choice of A. □ 

Because of the analogy between Equations ^ and (J7]), we define: 

Definition 8. The functions Kf,Hf of Q are the Gaussian and mean curvatures of 
the pair (m, s), i.e. of the polyhedral surface m with respect to the Gauss image s. They 
are associated to the faces of m. 

Obviously, mean and Gaussian curvatures are only defined for faces of nonvanishing 
area. They are attached to the pair (to, s) in an affine invariant way. There is a further 
obvious analogy between the smooth and the discrete cases: The Gauss curvature is 
the quotient of (infinitesimal) corresponding areas in the Gauss image and the original 
surface. 

4.3. Existence of principal curvatures. Similar to the smooth theory, we in- 
troduce principal curvatures K\, k,2 of a face as the zeros of the quadratic polynomial 
x 2 — 2Hx + K, where H, K are the mean and Gaussian curvatures. We shall see that 
in "most" cases that polynomial indeed factorizes, so principal curvatures exist. The 
precise statement is as follows: 

Proposition 9. Consider a polyhedral surface m with Gauss image s, and assume that 
for each face f G F mean and Gaussian curvatures Hf, Kf are defined. Regarding the 
existence of principal curvatures k\j and K2j, we have the following statements: 

1. For a quadrilateral f , k\ / = K2j <^=> rn(f), s(f) are related by a similarity. If 
this is not the case, Kij exist <^=^ the vertices of m{f) or of s{f) lie on the boundary 
of their convex hull. 

2. Suppose some linear combination of the n-gons m{f), s{f) is the boundary cycle 
of a strictly convex polygon. Then Kij exist, and k± / = K2j <^=^ rn{f) and s{f) are 
related by a similarity transform. 

3. Suppose f is a quadrilateral and the Gauss image s is inscribed in a strictly convex 
surface S. Then principal curvatures exist. They are equal if and only if m(f) and s(f) 
are related by a similarity transform. 

Proof. We consider the polynomial (p(x,y) := A(x ■ m(f) + y • s(f)) as in Prop. [6| 
The area of m{f) is nonzero, otherwise curvatures are not defined. Thus, </?(x, y) is 
proportional to ip(x,y) := x 2 — 2Hfxy + Kfy 2 ), and linear factors of ip correspond 
directly to linear factors of g(x) := tp(x, 1) = x 2 — 2Hx + K. So statements 1,2 follow 
directly from Prop. [6] As to the third statement, note that the vertices of an n-gon 
which lie in a planar section of X always are contained in the boundary of their convex 
hull, so we can apply 1. □ 
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Figure 2: Edge curvatures associated with a quadrilateral mo, . . . ,777.3 in a polyhedral 

surface m with Gauss image s. 



4.4. Edge curvatures. In a smooth surface, a tangent vector v £ T p M indicates 
a principal direction with principal curvature k, if and only if —dn{v) = kv. For a 
discrete surface m with combinatorics (V, E,F), a tangent vector is replaced by an edge 
(i,j) G i£. By construction, edges 7xi % 7xij are parallel to corresponding edges SjSj in the 
Gauss image mesh. We are therefore led to a curvature K e associated with the edge e, 
which is defined by 

(9) e = G E =>• Sj - Si = re^m* - mj) 

(see Figure [2]). For a quad-dominant mesh this interpretation of all edges as principal 
curvature directions is consistent with the fact that discrete surface normals adjacent to 
an edge are co- planar [16 . 

The newly constructed principal curvatures associated with edges are different from 
the previous ones, which are associated with faces. For a quadrilateral however, it is not 
difficult to relate the edge curvatures with the previously defined face curvatures: 

Proposition 10. Consider a polyhedral surface m with Gauss image s, and correspond- 
ing quadrilateral faces m(f) = (mo, • • • , "2,3), s(f) = (so, . . . , S3). Then mean and Gauss- 
ian curvatures of that face are computable from its four edge curvatures by 

/ lri N ^ ^01^23 — ^12^30 

v iU i = ; > 

«01 + ^23 - «12 - «30 

/ 1:J s K _ ^01^12^23^30 f 1 I 1 - —) 

«01 + ^23 - «12 - ^30 W2 K30 «01 ^23 ' 

Proof. We determine ctf such that x* := (mo V 777,2) H (mi V 777,3) = (1 — a/)wi + 
a/7773. Likewise we determine /?/ such that x* = (1 — Pf)rri2 + P/txiq. The condition 
So<i<4 K «,A;+i( m j+i — m i) = after some elementary manipulations leads to 

/-. v«23 + «30 KOl + ^12 7>^ /-i \ . 

if/ = (1 - a/) h a/ , Kf = (1 - a/)K23^30 + a/K iKi2- 

Hf = (1 -/?/) V Pf , Kf = (1 -P/)K30«01 +P/K12«23- 

Equating the two expressions for Hf and ify yields the result. □ 

Remark 2. Using the line congruence Lj = TXii V (mj + Sj) (cf. Section |2j, for each edge 
e = we define a center of curvature associated with an edge rxiirxij as the point 

c e = Li n Lj. The familiar concept of curvature as the inverse distance of the center of 
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curvature from the surface is reflected in the fact that the triangles OsiSj and c e mirrij 
are transformed into each other by a similarity transformation with factor l/ft e . 

5. Christoffel duality and discrete Koenigs nets 

We start with a general definition: 
Definition 11. Polyhedral surfaces m,s are Christoffel dual to each other, 

s = m* , 

if they are parallel, and their corresponding faces have vanishing mixed area (i.e., are or- 
thogonal with respect to the corresponding bilinear symmetric form). Polyhedral surfaces 
possessing Christoffel dual are called Koenigs nets. 

Duality is a symmetric relation, and obviously all meshes s dual to m form a lin- 
ear space. In the special case of quadrilateral faces, duality is recognized by a simple 
geometric condition: 

Theorem 12. (Dual quadrilaterals via mixed area) Two quadrilaterals P = (p±, 
Pi, P3, Pi) and Q = (q\, q2, q%, q^) with parallel corresponding edges, pi + \—pi \\ qi+i—qi, 
i £ Z (mod 4) are dual, i.e., 

A(P,Q) = 

if and only if their non- corresponding diagonals are parallel: 

(P1P3) II (<?2<74), (P2P4) || (qm)- 




Figure 3: Dual quadrilaterals. 



Proof. Denote the edges of the quadrilaterals P and Q as in Figure [3} For a quadri- 
lateral P with oriented edges a, 6, c, d we have 

A(P)= 1 -([a,b] + [c + d}), 

where [a, b] = det(a, b) is the area form in the plane. The area of the quadrilateral P+tQ 
is given by 

A(P + tQ) = -([a + ta*, b + tb*] + [c + tc*,d + td*}). 
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Identifying the linear terms in t and using the identity a + b + c + d = 0, we get 

4A(P,Q) = [a,b*) + [a*,b) + [c,d*) + [c*,d] 

= [a + b, b*} + [a*, a + b] + [c + d, d*} + [c* , c + d) 
= [a + b,b* -a* - d* + c*). 

Vanishing of the last expression is equivalent to the parallelism of the non-corresponding 
diagonals, (a + b) \\ (b* + c*). □ 

Theorem [12] shows that for quadrilateral surfaces our definition of Koenigs nets is 
equivalent to the one originally suggested in 0, [6] . For geometric properties of Koenigs 
nets we refer to these papers. It turns out that the class of Koenigs nets is invariant 
with respect to projective transformations. 



6. Polyhedral surfaces with constant curvature 

Let (m, s) be a polyhedral surface with its Gauss map as in Section [4] We define 
special classes of surfaces as in classical surface theory, the only difference being the fact 
that the Gauss map is not determined by the surface. The treatment is similar to the 
approach of relative differential geometry. 

We say that a pair (m, s) has constant mean (resp. Gaussian) curvature if the mean 
(resp. Gaussian) curvatures defined by ^ for all faces are equal. If the mean curvature 
vanishes identically, H = 0, then the pair (m, s) is called minimal. 

Although this definition refers to the Gauss map, the normalization of the length 
of s is irrelevant, and the notion of constant curvature nets is well defined for discrete 
surfaces equipped with line congruences. 

Theorem 13. A pair (m, s) is minimal if and only if m is a discrete Koenigs net and s 
is its Christoffel dual s = m* . 

Proof. We have the equivalence H = <^=^ A(m, s) = <^=^ s = m* . □ 

This result is analogous to the classical theorem of Christoffel [8] in the theory of 
smooth minimal surfaces. Figure [4] presents an example of a discrete minimal surface m 
constructed as the Christoffel dual of its Gauss image s, which is a discrete Koenigs net. 

The statement about surfaces with nonvanishing constant mean curvature resembles 
the corresponding facts of the classical theory. 
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Theorem 14. A pair (m,s) has constant mean curvature Hq if and only if m is a 
discrete Koenigs net and its parallel m l / H ° is the Christoffel dual of m: 

1 

m = m + — — s. 

Ho 

The mean curvature of this parallel surface (m + H ~ 1 s,—s) (with the reversed Gauss 
map) is also constant and equal to Ho. The mid-surface m + (2Ho)~ 1 s has constant 
positive Gaussian curvature Kq = 4Hq with respect to the same Gauss map s. 

Proof. We have the equivalence 

1 \ _ . 1 



A(m, s) = —HoA(m) Aim, m + -^-s) = m* = m + -r— 

V H o ) Ho 

For the Gaussian curvature of the mid-surface we get 



s. 



-2 

A i m +m^ A (rn) + ^-A(m,s) + (^A(s) 

□ 

It turns out that all surfaces parallel to a surface with constant curvature have 
remarkable curvature properties, in complete analogy to the classical surface theory. In 
particular they are linear Weingarten (For circular surfaces this was shown in |20j). 

Theorem 15. Let (m,s) be a polyhedral surface with constant mean curvature and its 
Gauss map. Consider the family of parallel surfaces m* = m + ts. Then for any t the 
pair (m ,s) is linear Weingarten, i.e., its mean and Gaussian curvatures Ht and Kt 
satisfy a linear relation 

(12) aH t + (3K t = l 

with constant coefficients a, (3. 

Proof. Denote by H and K the curvatures of the basic surface (m, s) with constant 
mean curvature. Let us compute the curvatures Ht and Kt of the parallel surface (m + 
ts, s). We have 

A(m + (t + S)s) _ 1 - 2H(t + 5)+ K(t + 5) 2 
A(m + ts) ~ 1 - 2Ht + Kt 2 

H-Kt 2 K 

l-2Ht + Kt 2 1 - 2Ht + Kt 2 
The last identity treats m + (t + S)s as a parallel surface of m + ts. Thus 



1 - 26 . "„ „., + S 2 . . = 1 - 2H,S + K,S\ 



m -. H -. Kt k, 



l-2Ht + Kt 21 1 - 2Ht + Kt 2 ' 

Note that H is independent of the face, whereas K is varying. Therefore, with the above 



values for Ht and Kt, relation (12) is equivalent to ]_-2Ht = = ^' wruc h implies 

1 t 9 

a=--2t, (3 = --t 2 . 



□ 
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Figure 5: A Koebe polyhedron s. The tangent cone from each vertex Sj touches S 2 along a 
circle Cj. These circles form a packing, touching each other in the points where the edges touch 
S 2 . It follows that the edge lengths are related to the opening angles oji of said cones: We have 

1 1 Si — Sj\\ — COt UJi + COt LOj . 



We see that any discrete Koenigs net m can be extended to a minimal or to a constant 
mean curvature net by an appropriate choice of the Gauss map s. Indeed, 

(m, s) is minimal for s = m*; 

(m, s) has constant mean curvature for s = m* — m. 

However, s defined in such generality can lead us too far away from the smooth theory. 
It is natural to look for additional requirements which bring it closer to the Gauss map 
of a surface. These are exactly three cases of special Gauss images of Section [2] 

Cases with canonical Gauss image. For a polyhedral surface m which has a face 
offset m! at distance d > (i.e., m is a conical mesh) the Gauss image s = (m 1 — m)/d 
is uniquely defined even without knowledge of m', provided consistent orientation is 
possible. This is because s is tangentially circumscribed to S 2 and there is only one way 
we can parallel translate the faces of m such that they are in oriented contact with S 2 . 
The same is true if m has an edge offset, because an n-tuple of edges emanating from a 
vertex (n > 3) can be parallel translated in only one way so as to touch S 2 . 

It follows that for both cases a canonical Gauss image and canonical curvatures are 
defined. In case of an edge offset much more is known about the geometry of s. E.g. 
we can express the edge length of s in terms of data read off from m (see Figure [5]) . 
The edges emanating from a vertex Sj are contained in Sj's tangent cone, which has 
some opening angle Uj. By parallelity of edges we can determine Uj from the mesh m 
alone. The ratio between edge length in the mesh and edge length in the Gauss image 
determines the curvature: Kij = ±(cot u>i + cot Uj)/\\mi — m,j\\ (we skip discussion of the 
sign) . 

7. Curvature of principal contact element nets. Circular minimal and cmc 

surfaces 

In this section we are dealing with the case when the Gauss image s lies in the two-sphere 
S 2 , i.e., is of unit length, ||s|| = 1. Our main example is the case of quadrilateral surfaces 
with regular combinatorics, called Q-nets. In this case a polyhedral surface m with its 
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parallel Gauss map s is described by a map 

(m,s) : Z 2 -> R 3 x S 2 . 

It can be canonically identified with a contact element net 

(m,V) : Z 2 — ► {contact elements in M 3 }, 

where is the oriented plane orthogonal s. We will call the pair (m, s) also a contact 
element net. Recall that according to [5] a contact element net is called principal if 
neighboring contact elements (m, T 3 ) share a common touching sphere. This condition 
is equivalent to the existence of focal points for all elementary edges (n, n') of the lattice 
Z 2 3 n,n', which are solutions to 

(m + ts)(n) = (m + ts)(n') 

for some 

Theorem 16. Let m : Z 2 —* M 3 6e a Q-net with a parallel unit Gauss map s : Z 2 — ► S" 2 . 
TTien m is circular, and (m, s) is a principal contact element net. Conversely, for a 
principal contact element net (m, s), the net m is circular and s is a parallel Gauss map 
of m. 

Proof. The circularity of m follows from the simple fact that any quadrilateral with 
edges parallel to the edges of a circular quadrilateral is also circular. Consider an ele- 
mentary cube built by two parallel quadrilaterals of the nets m and m + s. All the side 
faces of this cube are trapezoids, which implies that the contact element net (m, s) is 
principal. □ 




Figure 6: Parallel Q-nets m and m + s with the unit Gauss map s. All the nets are circular. 
The pair (m, s) constitutes a principal contact clement net. 



The mean and the Gauss curvatures of the principal contact element nets (m, s) are 
defined by formulas 

Proposition [9] obviously implies: 

Corollary 17. For a circular quad mesh m, principal curvatures exist w.r.t. any Gauss 
image s inscribed in S 2 . 

Recall also that circular Koenigs nets are identified in [7], |6] as the discrete isothermic 
surfaces defined originally in [3] as circular nets with factorizable cross-ratios. 
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Both minimal and constant mean curvature principal contact element nets are de- 
fined as in Section |6l It is remarkable that the classes of circular minimal and cmc sur- 
faces which are obtained via our definition of mean curvature turn out to be equivalent to 
the corresponding classes originally defined as special isothermic surfaces characterized 
by their Christoffel transformations. Since circular Koenigs nets are isothermic nets, 



from Theorem 13 we recover the original definition of discrete minimal surfaces from [3 . 

Corollary 18. A principal contact element net (m, s) : Z 2 -> I 3 x 5 2 is minimal if 
and only if the net s : Z 2 — > S 2 is isothermic and m = s* is its Christoffel dual. 

Similarly, Theorem [14] in the circular case implies that the discrete surfaces with 
constant mean curvature of |10|, |4] fit into our framework. 

Corollary 19. A principal contact element net (m,s) : Z 2 — ► M 3 x S 2 has constant 
mean curvature Hq ^ if and only if the circular net m is isothermic and there exists its 
dual discrete isothermic surface m* : Z 2 — > M 3 at constant distance \m — m*\ = The 
unit Gauss map s which determines the principal contact element net (m, s) is given by 

(13) s = H {m* -m). 

The principal contact element net of the parallel surface (m + j^, —s) also has constant 
mean curvature Hq. The mid-surface { m ~\"!^^ s ) has constant Gaussian curvature 4:Hq. 

Proof. Only the "if" part of the claim may require some additional consideration. If 
the discrete isothermic surfaces m and m* are at constant distance 1/Hq, then the map s 
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defined by (13) maps into S 2 and is thus circular. Again, as in the proof of Theorem 
this implies that the contact element net (m, s) is principal. Its mean curvature is given 
by 

A(m,s) A(m, Ho(m* — m)) 
A(m,m) A{m,m) '' 

□ 

7.1. Minimal s-isothermic surfaces. We now turn our attention to the discrete 
minimal surfaces m of [T]. which arise by a Christoffel duality from a polyhedron s which 
is midscribed to a sphere (a Koebe polyhedron). As Koebe polyhedra are up to Mobius 
transformations determined by their combinatorics, a passage to the limit allows us to 
determine in this way the shape of smooth minimal surface from the combinatorics of 
the Gauss image of its network of principal curvature lines. 

The Christoffel duality construction of pQ is applied to each face of s separately. We 
consider a polygon P = (po, . . . ,p n -i) with n even and incircle of radius p. We introduce 
the points qi where the edge pi~ipi touches the incircle and identify the plane of P with 
the complex numbers. In the notation of Figure [7] the passage to the dual polygon P* is 
effected by changing the vectors a« = q2i -z,bi = <?2i+i - z, a\ = p 2 i - $2i+i, b\ =p i — q i . 
Apart from multiplication with the factor ±/3 2 , the corresponding vectors which define 
P* are given by 

(i4) a * = (-iy/aj, b * = -(-iy/b-, a'* = (-iy/¥ p b '* = -(-iy/¥ r 

The sign in the factor ±/9 2 depends on a certain labeling of vertices. The consistency 
of this construction and the passage to a branched covering in the case of odd n is 




discussed in pQ. For us it is important that both P and P* occur as concatenation of 
quadrilaterals: 

(15) Pj = (j>j-iqjPjqj) for j = 0, ... ,n — 1 => P = Pi . . . © P„_i, 

and the same for the starred (dual) entities. The main result is the following: 

Theorem 20. A discrete s-isothermic minimal surface m according to pQ (Christoffel 
dual of a Koebe polyhedron s ) has vanishing mean curvature. Every face f has principal 
curvatures «i/,fv2/ = ~ K i,f- 

Proof. We first show that for all j, A{Pj,P*) = 0. This can be derived from pQ where 
it is shown that Pj and P? are dual quads in the sense of discrete isothermic surfaces 
[5] . Discrete isothermic surfaces are circular Koenigs nets jB], i.e., the quadrilaterals Pj 
and Pj are Christoffel dual in the sense of Definition [TTj 

We can see this also in an elementary way which tor p = 1 is illustrated by Figure 
[7J The angle a, = <(qi,z,pi) occurs also in the isosceles triangle q*z*q* +l , so non- 



corresponding diagonals in Pi,P* are parallel. By Theorem 12 A(Pi,P*) = 0. 

Lemma Snow implies that A(P, P*) = ^2A(Pj,P*) = 0. Thus all faces of m (i.e., 
the P*'s of the previous discussion) have vanishing mixed area with respect to s. As the 
faces of s are strictly convex, Prop. [9] shows that principal curvatures exist. □ 

7.2. Discrete surfaces of rotational symmetry. It is not difficult to impose the 
condition of constant mean or Gaussian curvature on discrete surfaces with rotational 
symmetry. In the following we briefly discuss this interesting class of examples. 

We first consider quadrilateral meshes with regular grid combinatorics generated by 
iteratively applying a rotation about the z axis to a meridian polygon contained in the 
xz plane. Such surfaces have e.g. been considered by |12] . 

The vertices of the meridian polygon are assumed to have coordinates (rj,0, h{), 
where % is the running index. The Gauss image of this polyhedral surface shall be gen- 
erated in the same way, from the polygon with vertices (r*,0,h*). Note that parallelity 
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Figure 8: Left: A polyhedral surface m which is a minimal surface w.r.t. to the Gaussian 
image s. Right: A polyhedral surface m! of constant Gaussian curvature w.r.t. the Gauss image 
s' (discrete pseudosphere) . 



implies 

rj+i ~n _ r* +1 



^ ^ hi + \ — hi h* +1 — h* 

Figure [8] illustrates such surfaces. All faces being trapezoids, it is elementary to compute 
mean and Gaussian curvatures of the faces bounded by the i-th and (i + l)-st 

parallel. It turns out that the angle of rotation is irrelevant for the curvatures: 



:y* . y r y» , iy* ~ -:y 1 ^ 

(17) H^= l \ 1+1 K« = ^±i 



*2 *2 

i 



1 2 > 1 1 



The principal curvatures associated with these faces have the values 
(18) k) = , n\ = . 

n+i + n r i+1 - n 

The interesting fact about these formulae is that the coordinates hi do not occur in 
them. Any functional relation involving the curvatures, and especially a constant value 
of any of the curvatures, leads to a difference equation for (ri)jgz. For example, given an 
arbitrary Gauss image (r*,0,/i*) and the mean curvature function defined on the 
faces (which are canonically associated with the edges of the meridian curve) the values 



rj of the surface are determined by the difference equation (17) an an initial value tq. 



Further the values hi follow from the parallelity condition (16). 

A meridian curve of a smooth surface of revolution does not intersect the rotation 
axis, and the Gauss map is spherical. Discrete analogues of such surfaces with a Gauss 
map s E S 2 and prescribed curvature are determined by the values (h*)i & z lying in the 
interval (—1, 1), and an initial value tq. The values r* = (1 — h* 2 ) 1 / 2 should be chosen 
positive. 

Remark 3. The generation of a surface m and its Gauss image s by applying k-th powers 
of the same rotation to a meridian polygon (assuming axes of m and s are aligned) is 
a special case of applying a sequence of affine mappings, each of which leaves the axis 



A CURVATURE THEORY FOR DISCRETE SURFACES 



17 



fixed. It is easy to see that Equations (17) and (18) are true also in this more general 
case. 



Remark 4. While the formula for K2 given by ( 18 ) is the usual definition of curvature for 
a planar curve, the formula for k± can be interpreted as Meusnier's theorem. This is seen 
as follows: The curvature of the z-th parallel circle is given some average value of 1/r 
(in this case, the harmonic mean of 1/rj and l/r r+ i). The sine of the angle a enclosed 
by the parallel's plane and the face under consideration is given by an average value of 
r* (this time, an arithmetic mean). By Meusnier, the normal curvature "sin a - ^" of the 



parallel equals the principal curvature kj, in accordance with (18). 



Example 1. The mean curvature of faces given by (17) vanishes if and only if ri + \ : = 
r* : r* +1 . This condition is converted into the first order difference equation 

(19) A In r; = -A lnr* (i e Z), 

where A is the forward difference operator. It is not difficult to see that the corresponding 
differential equation (lnr)' = — (lnr*)' is fulfilled by the catenoid: With the meridian 
(i,coshi) and the unit normal vector (— tanhi, 1/coshi) we have r{t) = cosht and 



r*(t) = 1/coshi. We therefore like to denote discrete surfaces fulfilling (19) discrete 
catenoids (see Figure [8] left). 

Example 2. A discrete surface of constant Gaussian curvature K obeys the difference 
equation KA(rf) = A(r* 2 ). Fi gurepl right illustrates a solutions. 



7.3. Discrete surfaces of rotational symmetry with constant mean cur- 
vature and elliptic billiards. There exists a nice geometric construction of discrete 
surfaces of rotational symmetry with constant mean curvature, which we obtained jointly 
with Tim Hoffmann. This is a discrete version of the classical Delaunay rolling ellipse 
construction for surfaces of revolution with constant mean curvature (Delaunay surfaces). 

Play an extrinsic billiard around an ellipse E. A trajectory is a polygonal curve 
P\,P2, ■ ■ ■ such that the intervals [Pi, Pi+i] touch the ellipse E and consecutive triples 
of vertices Pj_i, Pj, Pj+i are not collinear (see Figure [9J. Let us connect the vertices 
Pi to the focal point B, and roll the trajectory Pi, P2, . . . to a straight line £, mapping 
the triangles PPjPj+i of Figure [9] isometrically to the triangles PjPjPj+i of Figure 10 



We use the same notations for the vertices of the billiard trajectory and their images 
on the straight line, and the points Bi are chosen in the same half-plane of I. Thus 
we have constructed a polygonal curve B\,B^... . Applying the same construction to 
the second focal point A we obtain another polygonal curve A\,A%, . . ., chosen to lie in 
another half-plane of t. 

Let us consider discrete surfaces m and fh with rotational symmetry axis I generated 
by the meridian polygons constructed above: m, = Bi, fhi = Ai. They are circular 
surfaces which one can provide with the same Gauss map Sj := rrii — fhi. 

Theorem 21. Let Pi, P2, . . . be a trajectory of an extrinsic elliptic billiard with the focal 
points A, B. Let m, fh be the circular surfaces with rotational symmetry generated by the 
discrete rolling ellipse construction in Figures\§, 10- mi = Bi, fhi = Ai. Both surfaces 



(m, s) and (m, —s) with the Gauss map s = m — m have constant mean curvature H, 
where 1/H = \ A\B\ equals twice the major axis of the ellipse (see Figure^. 




Proof. The sum of the distances from a point of an ellipse to the focal points is inde- 
pendent of the point, i.e., 

I ■= \ABil 

is independent of i. Due to the equal angle lemma of Figure [TT] we have equal angles 
:= APxP-iAx = ZBP 2 P 3 and 7 := AP X P 2 B = ZP 3 P 2 A 2 in Figure [9j Thus P 2 in Fig- 
ure 10 is the intersection point of the straight lines {AiB 2 ) Pi {B\A 2 ). Similar triangles 
/\P 2 A\A 2 ~ AP 2 B 2 Bi imply parallel edges (AiA 2 ) \\ [B\B 2 ). This yields the propor- 
tionality rj/r, + i = r^ +1 /r- for the distances r to the axis i. For the mean curvature of 
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Figure 11: The angles between the tangent directions and the directions to the focal points of 
an ellipse are equal. 



the surface m with the Gauss image s = m — m we obtain from (17): 

lr i (r' i -r i )-r i+1 (r , i+1 -r i+ i) 1 



H 



•i+i 



The surface m is the parallel cmc surface of Corollary 19 



□ 



If the vertices of the trajectory Pi, P2, . . . lie on an ellipse E' confocal with E, then 
it is a classical reflection billiard in the ellipse E 1 (see for example |23j ). The sum 



d 



\APA + \BPA 



is independent of*. The quadrilaterals AiAi + iBi + \Bi in Figure 10 have equal diagonals, 
i.e., are trapezoids. The product of the lengths of their parallel edges is independent of 
i: 



(20) 



A4 A 



\BiB. 



i-Oi+l 



As we have shown in the proof of Theorem 21 is another product independent of 



i. An elementary computation gives the same result for the cross-ratios of a faces of the 
discrete surfaces m and fh: 



1 I A- A 



\B,:B. 



i-Dj+ll 



sin 2 a 



where 2a is the rotation symmetry angle of the surface. We see that q is the same for 
all faces of the surfaces m and fh. 

We have derived the main result of jllj . 

Corollary 22. Let Pi,P%,... be a trajectory of a classical reflection elliptic billiard, 
and m, fh be the discrete surfaces with rotational symmetry generated by the discrete 
rolling ellipse construction as in Theorem 21 Both these surfaces have constant mean 
curvature and constant cross-ratio of their faces. 

The discrete rolling construction applied to hyperbolic billiards also generates dis- 
crete cmc surfaces with rotational symmetry. 
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8. Concluding remarks 

We would like to mention some topics of future research. We have treated curvatures 
of faces and of edges. It would be desirable to extend the developed theory to define 
curvature also at vertices. A large area of research is to extend the present theory to 
the semidiscrete surfaces which have recently found attention in the geometry processing 
community, and where initial results have already been obtained. 
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